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Abstract 

We study quantum chromodynamics from the viewpoint of untruncated Dyson-Schwinger equa- 
tions turned to an ordinary differential equation for the gluon anomalous dimension. This non- 
linear equation is parameterized by a function P{x) which is unknown beyond perturbation theory. 
Still, very mild assumptions on P{x) lead to stringent restrictions for possible solutions to Dyson- 
Schwinger equations. 

We establish that the theory must have asymptotic freedom beyond perturbation theory and 
also investigate the low energy regime and the possibility for a mass gap in the asymptotically free 
theory. 
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1 Introduction 



We study non-perturbative aspects of quantum chromodynamics (QCD). We do so by investigating 
Dyson-Schwinger equations. Instead of solving a truncated version of these a priori very intricate 
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equations [|22ll . we use recent insight into the mathematical structure of quantum field theory to gain 
insight into the possible structure of solutions. This approach has been successfully applied to quantum 
electrodynamics in [5| and is here extended to QCD. 



1.1 The method 

We follow the methods employed in our work on QED [5], adopted to the study of QCD in the back- 
ground field approach as developed by Abbott [HI 121 • Let us first reconsider the situation for quantum 
electrodynamics. There, thanks to the Ward identity, it suffices to consider the anomalous dimension 
7i(x) of the photon which is essentially the /9-function, P{x) = X7i(x). 

This anomalous dimension is obtained from the photon's self-energy, a two-point function which is 
determined non-perturbatively by the knowledge of a single Lorentz scalar function 

oo 

G{x,L) = l-J2lk{x)L' , 

k=l 

with X the fine structure constant and L = In (— g^/ /i^), where the inverse photon-propagator is {q^g^y — 
q^qy)G{x,L). 

Combining this with the combinatorial Dyson-Schwinger equations and using an expansion into suit- 
able integral kernels which parametrize the corresponding integral equation, the Dyson-Schwinger 
equation combine with the renormalization group equation to give 

i) a recursion for the 7^: 

1k{x) = -77i(x)(l - xd^)'yk-i{x), k>2 , (1) 
k 

ii) a differential equation for 71 (x): 

7i(x)(l - x9,)7i(x) + 71 (x) - P{x) = . (2) 

Here, i) comes from the renormalization group, ii) from the DSE, and P(x) is a suitably constructed 
series over residues. The appearance of the operator (1 ± xd^) is typical for a gauge theory. 

We will now address a non-abelian gauge theory, resulting in a similar set-up (in particular, i) and ii) 
remain form-invariant) once we learned to make effective use of quantum gauge invariance to reduce 
again to a single ODE. Obviously, however, P(x) will be a different function, in particular, it will 
change sign. We study the consequences of this fact with minimal knowledge on the behavior of P{x). 
Still, as in QED, we will see that we can learn quite a bit regarding the non-perturbative sector of QCD. 

To proceed, we turn to the background field method. 
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1.2 Background field method 



We first have to consider the set of vertices and propagators in the background field gauge. They define 
a set TZ given as follows: 




cAc 



dAQQ 




\ / 



ccQA 






AQAQ 



Here, names attached to the vertices and propagators are short hand memos of the corresponding mono- 
mials in the Lagrangian. 

The field content is Q for the internal quantized gauge field, A for an external background gauge 
field, c, c for the anti-ghost and ghost field. Coupling of Fermionic matter will not change the ensuing 
discussion in any way and is omitted for convenience. See Abbott [j2J for details. 

With the set Tl comes an accompanying set of IPI Feynman graphs naturally labeled by elements in 
this set according to their type and number of external legs. 



We consider in particular Green functions for such graphs and adopt the results of 11181 . see also [[191 , 
which read as expected Vr G 7^ 



k>\ \-l\=k 



sym(7) 



B 



7;r 



ni;g7[01 ^ 
fl ^ 



£7(11 



with 



(3) 
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where maxf(r) is the number of maximal forests of F, |/i|v is the number of distinct graphs obtainable 
by permuting edges of h, bij(7, h, T) is the number of bijections of external edges of h with an insertion 
place in 7 such that the result is T, and finally (7|/i) is the number of insertion places for /i in 7 [|T8l . 
Sre<r> indicates a sum over the linear span < F > of generators of H. 

Next, we divide by the ideal / which implements the Slavnov-Taylor identities which here is generated 
order in order in (7^ by 

j^AAQQ ^ ^AQAQ ^ j^AQQQ ^ j^QQQQ ^ j^dQQQ 

On H/I we then get two independent Green functions which need renormalization 
corresponding to a mere two-element set 

TZh/i = {AA, cc} . 
Also, we then find an combinatorial invariant charge uniquely defined as 

C = I/VX^ = I/VXQQ , 
so that, as in QED, the /3-function is just half the negative anomalous dimension of the gauge field. 

Note that the addition of massless fermions would just add in H/ 1 an element tjjtjj for the fermion 
self-energy but would not change the ideal or the invariant charge. 

The system of combinatorial Dyson Schwinger equations is then 

k 
k 

Note that this determines in terms of itself, while X'^'^ is a function of itself and X^'^. We write 
X"" = I — X]fc>i[fl'^]''Cfc' with cl G H/I the generators of a sub Hopf algebra [Ull HI given by all 
graphs which contribute to an amplitude r at a chosen order k. Note that the resolution of a Green 
functions into images of Hochschild one-cocycles B"^ is the mathematical equivalent of a resolution of 
all overlapping divergences into non-overlapping integral kernels. That this is possible in a non-abelian 
theory was realized early by Baker and Lee [7 ]. 

See [[181 [191 HH |26l for explicit examples how Hochschild cohomology and Hopf algebras relate. 

There hence is a quotient Hopf algebra Haa spanned by generators c^'^ of X^^. Similarly, going 
temporarily to the quotient Hopf algebra Hcc defined by X^'^ = I, c^^ = 0, we find that this is a 
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cocommutative Hopf algebra (which is obvious from setting C = T) and hence we find a factorization 
of groups into an abelian subgroup Spec(ifcc) and a normal subgroup Spec{HAA), 

Spec{H/I) = Spec{HAA) x Spec(ifgc) , 
corresponding to a short exact sequence which splits 

I ^ Spec{HAA) Spec(i///) ^ Spec(i/g,) ^ I . 

It is this factorization which allows us to compute the /3-function of QCD by an ODE for a single 
equation below. Note that the situation is similar to QED: there, the Ward identity allows for a similar 
semi-direct product structure between photon amplitudes and Fermionic matter. Gauge invariance then 
allows to eliminate all short-distance singularities in the abelian subgroup thanks to the work of Baker, 
Johnson and Willey, and one is left with the photon propagation as the only source of renormalization. 

Here, we can compute the /9-function from G^^{x, L), but would have to consider the full coupled 
system to determine G^^{x, L) (and G'^'^{x, L)), which we do not attempt here. 

Furthermore, note that the simplification at A; = 1 

This is typically for gauge theories and emphasizes that we are in a single equation situation with s = 1 
O. Terms always deliver a pure residue from their short-distance singularities, and these terms 

are intimately connected to fermion determinants. We will not pursue this connection any further here. 

The background field method is then suited to our approach as it allows us to compute the QCD beta 
function from a single ordinary differential equation. 

Indeed, a change of basis of primitives allows to reduce the application of Feynman rules to the study of 
one- variable Mellin transforms for the integral kernels for the above primitives B']f{l), and from there 
we can strictly follow the techniques of l|20l |2T1 |27ll to get to an single ordinary differential equation: 

+ 1i{xf - P{x) - x-ii{x)-i[{x) = . 

Here, P(x) is a suitable series over primitives. As always, we renormalize using a momentum scheme 
with subtractions at = /x^. That scheme is uniquely suited [i l4.l to our gauge-invariant non-perturbative 
approach. 



1.3 Qualitative properties of QCD 

From perturbative computations, asymptotic freedom is firmly established. We will establish it beyond 
perturbation theory below. Perturbatively, this is mainly a self-consistency statement: assuming that 
the QCD coupling constant is small, we approximate the theory by its loop expansion to a few orders. 
The resulting polynomial approximation to the beta function supports the claim of asymptotic freedom 
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in perturbation theory: < 0, < x < 1, hence at large momentum transfer lim_Q2_^oo tts(L) — > 0, 
L = ln(— Q^/ /i^). The coupling indeed becomes small in that limit. 

As usual, perturbation theory agrees well with observations: asymptotic freedom is a well-established 
experimental fact. 

Much more intricate is the study 



This is beyond the reach of perturbation theory. Nevertheless, different approaches point out that in 
that limit, the gluon propagator might turn to a constant, confirming an old suggestion of Cornwall that 
the free gluon develops in the interacting theory a momentum-dependent mass which vanishes at high 
energies, but turns to a non- vanishing constant in the limit —Q^ 0+. We study this behavior from 
our viewpoint in section [3l 

Results to this effect were already obtained by 

i) Lattice computations 

ii) numerical study of Dyson Schwinger equations truncated in a gauge invariant way |fT3l : 

iii) in the Gribov-Zwanziger formalism [,17J . 

Below, we reconsider the problem from a study of the possible structure of solutions of Dyson Schwinger 
equations. We want to establish asymptotic freedom beyond perturbation theory, and want to discuss 
to what extent a solution which exhibits asymptotic freedom can also exhibit a mass gap. 

Again, as in the case of QED, we find the most interesting solution to be a separatrix. Li the case of 
QCD, that separatrix is the only solution which has asymptotic freedom. 



2 Results 



In QCD, the Dyson-Schwinger equation for 71 is 




(4) 



We will assume that the primitive skeleton function satisfies the following assumptions: 



H3: 



HI: 



H2: 



P is a twice differentiable function on R+, with P(0) = 0, P'(0) < and P"(0) < 0. 
There exist x* such that P(x) > —j and P"{x) < (i.e. P is concave) on [0, x*]. 
The function P{x) satisfies P{x) < for all x > 0. 
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As in [|5||, we avoid the singularities of (HI) at 71 = and x = by specifying with an initial condition 
at X*, namely, 

dji{x) , 7i(a;) +7i(x)2 - P(x) 

= f{li{x),x) = , 7i(x)=7o. (5) 

da; 3^71(3;) 

This ensures that solutions of Q exist at least locally around x = x*. Though we will mainly look for 
solutions or (HI) with 71(0;*) < and a; > 0, we will occasionally comment on the 71 (x*) > case. 

In the QED case, we proved in ^ (see also Section [5] of the present paper) the existence of a unique 
value 7i(a;o) (the separatrix) separating solutions that exist globally for all x > xq from those that 
can only be continued up to a finite x^nux > Xq. As shown in Section [51 all solutions in QED can be 
continued as a; — > 0, differing there 'only' by a flat behavior ~ e~^. In QCD the situation is reversed: 
we will prove that all solutions starting at some appropriate x* can be continued as x 00, but that 
there is a unique value 7i(x*) that separates solutions that cannot be continued as x — > from those 
that can, which either satisfy 71 (0) = — 1 if 71 (x*) < 7i (x*) or 71 (0) = if 71 (x*) = 7i (x*). We call 
the solution 7i(x) that satisfies 7* (0) = the asymptotically free solution. 

We will also use more speculative hypotheses on P(x): 

SI: There exists p > such that P(x) = — cx^ + o(x^) as x 00. 
S2: There exist x^ > —pr^ such that P"{x) < for all x G [0, Xc] 
S3: There exist a (finite) interval [xc, Xd] with Xc > x* such that 



2z 



-d;z > 1 . 



S4: There exist finite xi and Xr such that P{xi) = P{xr) = —\ and P(x) > — | for all < x < x/ 
andx > Xr- The function P(x) does no? satisfy SI, S2 and S3, butrather Iim3;^oo(-P(a^) — -Poo) = 
limj.^00 xP'{x) = for some Poo > —\- 



Let us briefly comment on our logic here. The H1-H3 hypotheses are a bare minimum for the results 
we will present below. Within perturbation theory, we have 

P(x) = 7i(x) + 0{x^) = -/3ix - (32x'^ + 0{x^) (6) 

as X ^ 0, where —/Si and —f32 are the 1 and 2-loop coefficients of the f3 function, namely /3i = 9 and 
(32 = 64 for Uf = 6. As such, the hypotheses HI and H2 are reasonable. While it also follows from 
dll) that P(x) < at least for small values of x, extending this to all values of x is somewhat more 
speculative. 

As we will show below, if P(x) < — | on a 'sufficiently large interval', for instance if either SI, S2 or 
S3 hold, all solutions of ^ satisfy 71 (x) = — 1 for some x > x*. It then follows from H3 that they 
grow linearly as x — 00 if 

f°° P(z) 
V{P) = - / -^dz < 00 , 

Jx* ^ 
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and faster than linearly if = c>o. Incidentally, we showed in Q (see also Section [5] of the present 
paper), that the finiteness/infiniteness of 'D{P) was intimately linked with the existence/non-existence 
as x — > oo of solutions of the analogous of ([5]) for QED. It is striking to see that the same criterion 
distinguishes between different type of behavior in QCD as welfl. 

Note that an anomalous dimension growing at least linearly as x ^ oo leads to a Landau pole for 
the running coupling, and hence a serious obstacle to studying the infrared behavior of the Gluon 
propagator. Despite that, we will show in Section [XT] that this pole (if present) can be removed using 
an unsubtracted dispersion relation, see e.g. Il24l . and the Gluon propagator can still be studied in the 
infrared limit. 



In contrast, if 71 (x) is finite as x — > 00, we avoid the Landau pole, and can study the Gluon propagator 
without using dispersion relations. Such constant asymptotics for 71 can only happen if P{x) tends to 
a constant as x ^ 00. This motivates (part of) the hypothesis S4. Under that hypothesis, we will show 
that there is only one solution that satisfies 



lim 7i(x) = = 7. 



00 



If Poo = 0, we call that solution the confinement solution 7j(x), and if Pqo > 0, we call it a strong 
confinement solution. On physical grounds, the asymptotically free and (strong) confinement solutions 
need to be the same. Unfortunately, for generic P(x) satisfying HI, H2 and S4, these two solutions are 
different. We did not succeed in finding a sufficient condition on P(x) that guarantees both solutions 
are the same. Despite that, a necessary condition is certainly that P(x) makes at least one (small) 
excursion below while avoiding the S1-S3 conditions, see also figure[I]below or Section l44l 

By standard folklore and heuristics |l3l, a nowhere vanishing /3-function, /5(x) < 0, Vx > 0, which 
we will indeed establish below under assumption HI above, implies a mass gap in QCD, and hence a 
confinement scenario following old ideas of Cornwall [[TSll . We will come back to that in Section [3] of 
this paper. 



Remark 2.1 If Poo > in hypothesis S4, P(x) has a further zero for a finite xi > 0, P(xi) = 0. In 
such a case, using the running coupling formulation of0!jf, one sees that some solutions spiral around 
the zero of P{x), themselves having infinitely many zeroes. The asymptotically free solution 7^(x) may 
or may not spiral around x = xi, depending on details of P. Should it be captured, we get a solution 
71 (x) which has a UV fix-point at zero and an infrared fixpoint at xi. This is the Banks-Zaks scenario 

m. 



In the remainder of this section, we are going to state our main results. The proofs and technical details 
are postponed to Section |4] of this paper. 

Our first main result gives a complete characterization of the behavior of solutions of (|5]) for x < x* . 
In particular, it establishes the uniqueness of the solution exhibiting asymptotic freedom. 

'though of course the primitive skeleton functions P are different in both cases 
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Theorem 2.2 Under the hypotheses HI and H2, there is a unique value 7i (x*) < such that the 
corresponding solution -^l^x) o/© exists for all x G [0,x*] and satisfies lim7j'(x) = 0. Additionally, 

that solution satisfies 

^™'-^ <7;W<P'(0).. (7) 

If'ji^x*) < 7j'(x*), corresponding solution satisfies lim7i(a;) = —1. If'y^{x*) < 7i(a;*) < 0, 

then there exist Xmm > such that the corresponding solution satisfies 7i(a;min) = 0. 



Note that since ^^^^^^""^ ^ = P'{0)x + 0{x'^) as x ^ 0, © shows that ^^{x) = P'{0)x + 0{x'^) as 
X ^0. 

As we will prove in Proposition |4]8]) of Section|4] below, all solutions of Theorem |2.2| can be continued 
as X ^ oo by only adding the H3 assumptions. In Proposition I4.6[ we will show that solutions that 
satisfy 7i(xmin) = for some < Xmin < x* can be continued in the first quadrant (becoming double- 
valued) by reverting to the so-called 'running coupling' formulation of © (see also [|5|). In particular 
these solutions will satisfy 71 (xq) > for some xq > Xmin- 

Our second main result concerns the asymptotic behavior as x — 00 of solutions that enter the first 
quadrant, or attain the value —1 somewhere. This last condition can be verified under additional as- 
sumptions on P such as SI, S2 or S3. 

Proposition 2.3 Assume P(x) satisfies H1-H3 and that one of the two following statements holds: 

1. —1 < 7i(x*) < and P(x) satisfies SI or S3, 

2. 7i(x*) < 7i(x*) and P{x) satisfies S2. 

Then there exists Xq > x* such that the corresponding solution 71 (x) satisfies 71 (xq) = — 1. 

To be able to state our asymptotic result as x — »• 00, we need first to introduce the slope function 
Sp (xo , x) . This function is given by 

s.(.o.x)^ (21(^+2 rz^i)^. 

Note that, if V{P) < 00, the slope function Sp(xo, x) goes to a finite value as x ^ 00 for any xq. 
We can now completely describe the asymptotic behavior as x — 00 of solutions of Q: 
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Theorem 2.4 Assume P{x) satisfies H1-H3. If there exists Xq > such that either 71 (xq) > or 
7i(a^o) < — 1 then 

X S'p(xo,x) < 7i(x) < xi Sp{xo,x) + — //■ 7i(xo) > , 

-X S'p(xo,x) < 7i(x) < -X ( S'p(xo,x) - — ) - 1 if 1i{xq) < -1 . 

Furthermore, ifV{P) < 00 and '~^i{xq) < —1 or 7i(xo) > 0, ?/zere exists s > 5mc/z 

7i(x) f -s < if ^i{xo) < -1 

iim = < 

X y s > if 7i(xo) > 

7/'7i(xo) < —1, convergence towards the limit is given by 



< C r -^dz . (8) 



X 



X 

IfV{P) < 00 and 7i(xo) > 0, then ^ also hold, with — s replaced by s. 

We want to stress here that the slope value s depends on the actual solution under consideration. Also, 
for solutions that eventually enter the first quadrant, the value of the slope s has no reason to be the 
same along the two branches of the solution (the one in the first quadrant, and the one in the fourth). 
Also, note that the result depends only on the assumption 71 (xq) = — 1. 

We now state the existence and uniqueness of confinement solutions under hypothesis S4. 

Theorem 2.5 Assume P{x) satisfies HI, H2 and S4. Then there exist a unique solution 7j(x) 0/(0) 
such that 

r CI , r 1 + Vl + 4P(x) 

hm 7i(x) = hm = 700 . 

a;— >oo x— >oo Z 

Furthermore, 7f (x) satisfies the usual trichotomy as x decreases: either 7i(0) = —1, or 7i(0) = 0, or 
^l{x) cannot be continued for x < x^^^ for some x^^^ > where 7i(xmin) = 0. Finally, if there exists 
a^max such that 'P(x) > 0/or all x > Xmax (hence P(x) is strictly increasing towards Poo), then 



v/r+4i^+l ^^ Vl+4P(X) + 1 

7: < < 



for all X > Xmax- 

As already noted above, we cannot show that 7i(x) = 7j(x) without additional hypotheses on P(x). 
We can however note that the two types of solutions are compatible: the confinement solution satisfies 
the integral equation 
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Figure 1: A 'fragile' solution: the solid bold curve is the solution combining asymptotic freedom 
and confinement for some artificial P(x) with Pqo = 0. The dashed curves are the nuUclines 71 + 
7^ — P{x) = 0. The solid curves 7j_|_ and 7j'_|_ correspond to confinement and asymptotically free 
solutions with P{x) perturbed so that the gap \xr — xi \ is slightly larger ('+' subscripts) or smaller ('—' 
subscripts). In both cases, the asymptotically free solution and the confinement solution do not match. 



whose r.h.s. converges to as x — > if P(a;)/7j(x) ^ 1 as x — 0. Characterizing the set of functions 
P(x) for which the confinement solution and the asymptotically free one are the same is a difficult 
problem, see e.g. figure [T] for an example with an 'artificial' P{x). A necessary condition is that 
P{x) makes an excursion below — | on (at least) one interval so that the nuUclines, i.e. the location in 
(x, 7i) -plane where 7i(x) = 0, 



±v/l+4P X 
l±{x) = — 



show a gap as in figure [T] However if the gap is 'too wide', we have 7j'(x) = 700 for some finite x, 
and 7i(xmin) = for some Xmin > 0, whereas if the gap is 'too small', 7i (x) intersects the nuUcline 
at some x > Xr, and hence cannot reach 700 as x — > 00, while 7i(x) intersects the nuUcline at some 
X < xz, and hence cannot reach as x ^ 0. 
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We conclude this section by explaining the terminology asymptotically free, confinement and strong 
confinement for our solutions. These come from the study of the inverse Gluon propagator, 

p-i(x, Q2) = Q^G{x, L) with L = ln(gV/"^) • 

In Section (|3.2I) below, we solve the RGE equation expressing the scale invariance of G{x,L). In 
particular, in Theorem 13 .21 we show that if 71 (a;) 700 as x ^ 00, then 

X 

where X(t, x) is the running coupling, i.e. the solution of 

dX{t,x) ^x(t,xhi(X(t,x)) with X(t = 0,x) = x. 
at 

The possible large \t\ behavior of X(t, x) are given by 

XooG"'°°^ as t ^ —00 if lim7i(a;)=7. 



00 ; 

x^oo 



X{t,x) ~ < 



Xoe~* as t ^ 00 if lim7i(x) = —1 , 

x—>0 



where Xoo and xq are some positive functions of x. Thus, in all cases where 71 (x) can be continued to 
X = 0, G(x, L) — > in the ultraviolet regime L ^ 00, and moreover 



P-\x,Q' 



) ^ 

+xo as L — s> 00 if lim7i(x) = — 1 , 



as L ^ 00 if Um2lM = p'fQ) 



x^O 



hence the terminology asymptotic freedom for the solution that satisfies 71 (x) = xP'(O) + O(x^) as 
X — s> 0. In the infrared regime, however, 

G(a;,L)^^e^-^ = ^f^V°° as L -00 , 

X X \fX^ J 

and hence we find 

^2 jf jjj^ 7i(x) = -1 

lim P'\x,Q')={ ' ™ ^ ^ , (10) 



L— >— 00 



00 if lim 7i(x) < —1 



a finite mass gap if 7c 
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3 The Gluon propagator, confinement and mass gaps in QCD 



3.1 Corrections from dispersion relations 

In QED, the coupling is weak at low energy or momentum transfer. In our previous work fS^, we 
could hence define boundary conditions at low energy, and studied the behavior of J2k IkL^ for L > 0. 
In particular, for 71 = 71 (a; (L)), a continuation to L < was never needed by the choice of our 
renormalization conditions. On the other hand, for large L ^ we could establish a separatrix, which 
possibly avoids a Landau pole at any finite positive L. We conjectured that this might be the solution 
chosen by Nature, and further detailed analysis of its properties awaits more analysis of the function 
P{x). Should it turn out that P{x) is such that the separatrix will not avoid a Landau pole {i.e. turns 
to infinity at finite L), we will have to turn to dispersion relations to understand the non-perturbative 
corrections coming with such a pole, as recognized by Shirkov and collaborators early on ifTTI . 

For QCD, we again fix a small coupling, but this time large momentum transfer, L for our 
boundary conditions. We are now interested in a continuation to L ^ 0, in particular we are interested 
in L ^ —00. Under very mild assumptions on P{x), and certainly by any experience from perturbative 
approximations of the theory, we expect the anomalous dimension 71 to go below the value —1 at some 
finite coupling x, and hence x{L) to turn to infinity at some finite negative L. Shifting that L to zero 
essentially defines the scale Aqcd, and we are interested for that shifted La to study the regime La < 0, 
in particular La —00. To consider such a limit based from an approach formulated for La > 0, we 
will use dispersion relations. Our approach is motivated again by Shirkov and collaborators work [|24ll . 

On general grounds, we know that x{Ljyj and G{x, La) can be treated by an unsubtracted dispersion 
relation [23J: 



The inverse propagator needs a subtracted dispersion relation, which leads back to an unsubtracted 
dispersion relation for G [|23l . 

The Dyson-Sch winger equations themselves are supposed to hold for the whole theory regardless of 
the sign of La. Similar, the renormalization group equations for the running coupling are supposed to 
hold. Our derivation which turned the Dyson-Schwinger equations into a ODE was valid for L > —La, 
hence remain valid, after shifting, for La > 0. 

Continuing to La < will generate non-perturbative corrections to 71(2;), and hence P{x), determined 
from the requirement that equations of motion, renormalization group flow and analyticity properties 
of field theory are what they are supposed to be. 

Any perturbative approximation is in accordance with these properties of field theory only up to the 
order considered. When we study solutions of Dyson-Schwinger equations, we demand accord with 
these properties as a guide to find the necessary non-perturbative corrections in the region La < 0. 
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We hence will start by first applying a dispersion relation to analyze x(La). That leads to a corrected 
7i,disp due to the fact that at La = we find that x = oo. Combining this with the assumption that the 
uncorrected 71 (x) is driven below -1 at finite x allows an easy estimate of the corrected 7i,disp and also 
allows for consistency with the direct analysis of G by dispersion methods. 

Let us now start with a study of x{L\). We start our considerations by boundary conditions such that 
G{x,L) = 1 at some very high momentum transfer ^ 0, L = ln((5^//i^), which determines a 
suitably small x in agreement with say deep inelastic scattering experiments IfTOl . 



dxix,L) _ , , 1 , 



We have 



dL u'yi{u 

and assume that the integral 

du < 00 



M7i(m) 



In particular, we assume that 71 (x) < —1 for some positive finite x, in accordance with our previous 
discussions and experimental evidence. We thus define 



,2 _ 2 , {/r^^^d"} 



Aqcd = + e 

Setting L\ = ln(Q^/AQQj-,) to absorb the dependence on x, jj?, we get 

1 



(La) 



du = LK. (11) 



This equation defines La{x) as well as the inverse function x(La). Since x(0) = +00, we cannot trust 
our solution for < A'^. To remedy this, we use our previous result that for large x, 71 (x) —sx, 
under our present assumptions. 



Following the conventions of Shirkov 11241 . we define a running coupling in accordance with the ex- 
pected analytic behavior of field theory using a dispersion relation: 

The pole at —Q"^ gives us back the uncorrected x(La) . But by assumption, there is a further pole in the 
complex (T-plane, located at La (00) = 0. To study the contribution from that pole, we first note that as 
X — s> 00, we have 

-au ~ / — - du 



i{LA) ulliu) Js^^La) sx(La) ' 

and hence x(La) ^ near La = 0. Feeding this relation into (fT2l) gives (see f2M ) 
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If we were to identify s with the one-loop coefficient of the /5-function, this would reproduce Shirkov's 
analysis for one-loop QCD, see [24], where Shirkov also notes that s seem not to vary much at low loop 
orders. Note that the correction to x{Ljyj goes to the finite value 1/s in the infrared limit La — oo. 

Now, 7i also obeys an unsubtracted dispersion relation. As 71 (x) is finite for all finite x and depends 
on L only through x, the dispersion integral will correct 7i(a;(L)) ^ —sx{La) by 

1 

7i,disp(a;(LA)) = 7i(x(La)) - ■Y~^ ■ 

Note now that the correction to 71 goes to —1 as La ^ —00. Using ([T]) and x'yi{x)dx = di, we get 

(-1)^ 

7fc,disp(a;(LA)) = 7fe(x(LA)) + 7fe,corr(-^^A) With 7fc,corr(-^^A) ^ ^, as La -00 . 
As such, the correction from the dispersion relation to the function G{x, La) satisfies 

Gcorr(a;, La) ^ S:j^^A = ^'^^ = rfi La ^-00, 

which gives an inverse propagator satisfying 

lim p-^(x,Q2) = -A^ 

This gives a finite and renormalization group invariant positivity-violating mass gap. This compares 
nicely with the results of Gracey et.al. ifTTl . 

Note that we rely completely on the assumption that physical quantities in massless field theory have 
neither poles nor branch cuts off the negative (in our conventions) real axis and the result that the 
solutions we get from our ODEs for the uncorrected x, 7^, G(x, La) are in accordance with these re- 
quirements but for the isolated pole at La = 0. 

If a future analysis of P{x) justifies the trust in field theory expressed in this section remains to be 
seen. We are content having identified a clean mechanism for the generation of a mass gap, based on 
the consequences of the Hopf algebra structure underlying local field theory, and assuming analyticity 
properties in accordance with the underlying axiomatic structure of local quantum fields. 

Let us add a few comments which put our results in context. We rely on the assumptions on the function 
P{x) outlined above. Note that the ODE for the anomalous dimension 71(2;) takes into account all the 
iterations of superficially divergent graphs into each other, hence all renormalon ambiguities met in the 
context of resummation of a perturbative series are taken care of. We hope this will make P{x), which 
has an asymptotic expansion related to a weighted skeleton expansion, amenable to more constructive 
methods of analysis in the future, in contrast to 71 (x) whose resummation as a perturbative series faces 
such ambiguities. 

On the other hand, we emphasize that our analysis avoids any truncation of Dyson-Schwinger equa- 
tions, and any assumptions made on the infrared behaviour and infrared powercounting of QCD ampli- 
tudes. It thus complements the approaches available so far in the literature, see for example 1. 16,1 . While 
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other methods based on Dyson-Schwinger equations attempt to solve them in some justifiable limit or 
truncation, we establish properties of the solution of the full equations depending on assumptions made 
for P(x). 



3.2 Solving The Renormalization Group Equation for G{x, L) 

Our goal in this section is to analyze the dressing function G{x, L) that modulates the free inverse 
propagator. On general grounds, this function solves the Renormalization Group equation 

^— 9l + x'-)i{x)dx — S7i(x) ^(^(x, L) = and G{x, 0) = 1 for x > , (13) 

where s = ±1 distinguishes between writing the propagator as 



with s = lors = — 1. Note that the difference of (fT3l) from the usual RGE equation 

>^ + /3(^)9, + 7(^?))g = 



is merely a matter of convention on the definition of (3, 71 and x. In particular, (1131) agrees with two-loop 
computation of the Gluon propagator. 

We will here make the following (minimal) hypothesis on 71 (x): 

Hypothesis 3.1 The function 7i(x) is a negative C^([0,oo), (— oo,0]) function whose only possible 
zero is at x = 0, where 71 (x) = —dx'^'^ + C(x*'^^) with qo >0 and d > 0. 

Note that the results of Section [21 give qo = I and d = — P'(0) for the asymptotically free solution 
7i(x), in accordance with perturbation theory. For the solutions satisfying 7i(0) = —1 however, we 
get go = and d = 1. 



We first note that one can attempt to solve (1131) by writing 

7x (x) 

ikyxju Willi ■ikyxj = 

k=l 

since one has (at least formally) that 



00 /X 

G(x,L) = l-s^7fc(x)L^ with 7,(x) = ^(x9,-s)7fc-i(x) k>2, (14) 



00 

Ol + X7i(x)9a; - S7i(x)j(^(x, L) = -S7i(x) + s ^ /i;7fc(x)L^~^ - 7i(x)(^x9^ - s^'jkix)L'' 

k=l 

(7fc(x) - ^(^^x - s)ik~,{x)\ kL'^ = 0. (15) 

k=2 ^ ^ 
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This approach naturally raises the (difficult) question of convergence of the series in (fT4l) and (fTSi) . 
We will use instead an alternative way of solving (fT3l) that avoids these convergence problems. In 
particular, while the series (fT4l) necessarily converge on a symmetric interval of the form {—Lq, Lq) for 
some (possibly infinite) Lq > 0, our approach will give a solution to (fT3l) that is defined on an interval 
of the form {—Lq, oo) for the same Lq. As such, the approach below shows that the limit L ^ oo of 
G{x, L) makes sense also if the series solution converges only on a finite interval. 

Our method is based on the fact that (fT3l) can be transformed into a linear transport equation by ap- 
propriately factorizing G{x, L) into one part that cancels the term involving no derivatives and one part 
that solves a genuine transport equation of the form 

{-dL + x^^{x)d^)H{x,L) = {] . (16) 

As such, it is important at first to consider the characteristics curve of (fT6l) . For each fixed x > 0, we 
first define X{t^x) as the solution of the running coupling equation 

dX{t,x) ^ x{t,x)-fi{X{t,x)) with X{t = 0,x) = x . (17) 

CJ.6 

Since 7i(x) is assumed to be C^, solutions of this equation exist at least locally around t = 0. For 
further reference, we denote by V{x) the maximal interval of existence of the solution of (fTTI) for a 
fixed X. 



Then for each fixed {x, L) E x R, we define the characteristic curve C{x, L) as 
C{x,L) = ^^{X{t,x),L-t) with teV{x)^ , 

X,L- [ -^^^ with XgR+1 . (18) 

Note that both above formulations of the characteristic curve C{x, L) are equivalent. The characteristics 
corresponding to different values of L are vertical translations of the same curve in the (x, L)-plane. 
By hypothesis 13. 1[ we get that the characteristics are all asymptotically vertical as X ^ 0, and, as a 
function of t, we have 

I e"*^* as t — i> oo if f/o = 
X{t,x)^{ 1 -f a- (19) 

I (ct) -90 as t ^ OO if go > 

The behavior of the characteristics as X ^ oo depends on the asymptotic behavior of 71 (x) as x ^ 00. 
In all cases, C(x, L) approaches (00, L — Loo(x)) as X ^ 00, where L^oix) is the possibly infinite 
quantity defined by 



-'00 



x 



dz 



X 2:71(2) 



This shows that the maximal interval of existence for solutions of (fTTI) is V{x) = (Loo(x), 00). The 
results of Section |2] show that Loo(x) is finite if there exists xq > where 71 (xq) = —1, and infinite 
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Figure 2: Generic shape of characteristics curves in (x, L)-plane. The left panel shows the case where 
Loo{x) is finite, the right panel when it is infinite. In the left panel, the characteristics cross the line L = 
if and only if they are above the (bold) curve L = Loo (a;), and all characteristics are asymptotically 
horizontal as x — > oo. In the right panel, all characteristics cross the line L = 0, as they link (0, —oo) 

to (CXD, oo). 



18 



otherwise. If Loo{x) > — oo, the characteristic C{x, L) intersects the line L = if and only if L > 
Loo{x). If —1 < 7i(x) < for all x > 0, then Loo{x) = — oo for all x > 0, and all characteristic 
curves cross the line L = 0. The generic shape of the characteristics curves is displayed in figure [21 

We can now give the solution of (fT3l) in accordance (for s = — 1, as expected) with Bogoliubov-Shirkov 
(Ea, App.IX, eq.(27)): 



Theorem 3.2 Assume 7i(x) satisfies Hypothesis \3.1\ Then the solution o/dT?]) is given by 
for all (x, L) such that Loo{x) < L < oo and x > 0. 



We want to stress here that the relation (|20|) hold only for pair of values of (x, L) satisfying L^{x) < 
L < oo. For pairs (x, L) with L < Loo(x), G(x, L) can only be determined from ^(x, 0) = 1 if one 
specifies 71 (x) for x < as well. 

Assuming X(t,x) to be analytic for all t G 2^(x), it is a straightforward computation (see Section 
(13.31) ) to show that (fT4|) is the Taylor series expansion of (|20|) at L = 0. As such, the series solution 
(fT4l) is expected not to converge for \L\ > Loo(x). However, the function provided by (|20l) is defined 
for unbounded positive L and do solve (fT3l) for all such L. If -Loo(x) > — 00, Theorem 13.21 thus gives 
^(x, L) for values of (x, L) for which the series formulation fails. 



proof of Theorem 13.21 We first set 



G{x,L) = x'H{x,L) 
for all X > 0. Substitution into (fT3l) gives 

^- 9l + X7i(x)9^jif(x, L) = and H{x,0) = x^"" if x>0. 

We then note that H{x, L) is constant along the characteristic curve C(x, L), for we have 

^(/f(X(t, x), L - t)) = -D2H{X{t, x),L-t) + DiH{X{t, x), L - 

c^L + X7i(x)(9^)f/'(x, L) 



= . 

x=X{t,x),L=L-t 

If the curve C(x, L) crosses the line L = in the (x, L)-plane, it does so at t = L, and we get 

H{x, L) = H{X{L, x), 0) = X(L, x)"" , 



(21) 



(22) 



which completes the proof. ■ 
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3.3 The series solution for G{x, L) 



Throughout this section, we consider 71 to be a fixed solution of (|5]) that exists for all x > 0. We first 
pick xo > and to G R- These values are arbitrary. We then introduce 



T{x) =to + 



(23) 



We then note that 



lim T{x) 

X— >oo 



T 



-00 if - 1 < 7i(x) < Vx > 



Too > — c>o if 3xo > with 71 (xq) 



(24) 



Since T'(x) = — 7^- < 0, x t— > T(x) is an invertible map from R+ to [Too, 00) • Moreover, we have 
T(0) = 00 and T{oo) = Too- The inverse map x(t) (the solution of t = T{x{t))) satisfies the 'running 
coupling equations' 



'^'''^'^ ~-Ut)+li{t?-Pm) and ^ = x(t)7i(t) 



(25) 



where 7i(t) = 7i(x(t)). Since T(oo) = Too, x{t) diverges as t — Too (and so does 71 (t) if 71 (x) = — 1 
for some x > 0). 

Consider now the series solution (fl?)) . We first introduce the functions Sk such that 

1 



k\ 



Sk{T{x)) 



for /c > 1. Substitution into (fT4l) gives 



S,{T{x)) = S',_,{T{x))={ ^^fc-iW 



t=T{x) 



Since s = ±1, we find from ^ that Si{t) = g|l = -sf^{x{t)-'). Using = 1, we find 



X* / d^ / 1 \ 



(26) 



f=T(x) 



for all k > 1. Since x(T(x)) = x, the r.h.s. of (|26l ) is equal to 1 when A; = 0, and so 



G(..L) = l + 5:-n.WL' = x-E ( ( )) 

fc=l fc=0 ^ \ J / 



(27) 



t=r(a:) 



We now fix < X < 00. Since x is finite, T(x) > Too, and since the above series is a Taylor series 
of x(t)^** at t = T(x), the series converges and takes the value x{t + T)^"* for small values of L if we 
assume z{t) to be analytic att = T(x). We thus find 



G{x,L) 



X 



x(T(x) + L) 
20 



(28) 



at least for sufficiently small L. This formula is the same as the one of Theorem |3.2t x{T{x) +t) solves 

'^''^'^^^^ = x(T(x) + t)7i(x(T(x) + t)) , (29) 

with initial condition x{T{x)) = x att = 0, and so x{T{x) + t) = X{t, x) by uniqueness of solutions 
of CI]) and (Ell). 



The question of convergence of (1271) for large L depends on the (in-)finiteness of Too. If T^o is finite, 
x{t) has a pole as t ^ Too, and (l27l) must diverge for 

L <Too- T{x) = / = Loo(a;) 



X 



as in Section |321 However, since the series (l27l) can only converge for L on symmetric intervals, it is 
also divergent for L > —L^{x) > 0. In that case, we have to revert to Section [X2l for the validity of 



4 Technical details for QCD 



In this section, we present the technical details proving the statements of Section [2] on QCD. In our 
analysis of ([5]), we use mainly two tools: integral representations of the solutions and nuU-clines. Our 
strategy is directly inspired from |i6J and our previous work [5] on QED. 

There are two types of integral equations one can write for dS]). The first one (see also ^) reads 



_ x(l + 7i(x-)) _ r P(z) 



liix) = - 'i^^-l-x -Y^dz. (30) 

For the second one, we let 71 (x) be a solution of ([5]) and 72(0;) be any function. Then for all xq, x G X, 
where X is the common interval where 71 and 72 are defined, we have 



71(3^) - 72(3;) = (^71 (a^o) - 72(3^0) j ^[71, 72](a;o, a;) + J R[-f2]iy) K[-fi,-f2]iy, x) dy , (31) 

where 

dx X72(x) ' 
^[71, 72] (3:0, 3;) = ( — ] exp( - [ P^^—dz] . 

The nuU-clines are defined as the locations in (x, 7i)-plane where solutions satisfy 7^ (x) = 0). These 
are given by the graph of the two functions 



21 



In particular, as -P(O) = by hypothesis HI, and P{x) > — ^ for x G [0, x*], the null-clines extend at 
least up to the line x = x* in the (x, 71) -plane. On the null-clines, the second derivative of 71 is given 
by 

ax ^ |x7^(x)| 

By hypotheses HI and H2, we have 

P'{x) = P'{0)+ [ P"{z)dz < P'{0) <0 yx e[0,x*] . (33) 
Jo 

Hence by (|32l) . solutions of ([5]) can have at most one local maximum in the interval x E [0, x*], and no 
local minimum. For further reference, we also note that by hypotheses HI and H2, 

P{x) = P'{0)x + (^j^ P"{z)dz^ dy < P'{0)x Wx e [0, x*] . (34) 



In particular, we have 

. -P(x) 



mm 



P'(0) = |P'(0)|>0, (35) 



xe[o,x*] X 

and since P{x*) > — | by hypothesis on x*, we have x* < —^prj^j < 00. 

4.1 Existence, uniqueness and properties of the asymptotic freedom solution 

We can now establish the existence, uniqueness and properties of the solution with asymptotic freedom. 



Theorem 4.1 Under the hypotheses HI and H2, there exists a unique value 7i (x*) such that the cor- 
responding solution 7j'(x) o/© exists for all x G [0, x*] and satisfies lim7j'(x) = 0. 

Proof. We first prove that the solution is unique. Namely, assume ab absurdum that 71 and 72 are two 
solutions of © on [0, x*] satisfying lim7j(x) = 0. Since 7i(x*) < and lim7j(x) = 0, we necessarily 

X— >0 x^O 

have 7i(x) > 7(f (x) for all x G [0, x*]. Since 72 is a solution, we can apply (|3T1) with ^[72] = 0, and 
using (|35l) . we get 

|7i(^) -72(a:)| > - ^,[x')\^^e^J f ^^dz) . (36) 

X \Jx iti^r J 

Since 7^^ = P'{Q)z + 0{z'^) as z ^ 0, the r.h.s. of (l36l) diverges as x ^ 0, which contradicts 
lim7j(x) = 0. This shows that there can be at most one solution of ® satisfying lim7i(x) = 0. 

IT— >o x— >o 
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To prove the existence of this solution, we define the two sets 



Ii = e]-f^{x^),0[s.t. 3s^in e]0,x*[with7i(xmij = 0} , 

h = G]7+(x*),0[s.t. 3xi G]0,x''[with7i(xi) = 7+(xi)} . 

We will prove in Propositions 14.21 and l4.4l below that these sets are non-empty. Continuity of solutions 
w.r.t. initial conditions imply that they are open, while (|3TI) with R['j2] = shows that solutions are 
ordered, which imply that each Ij is a single interval. Now since Ii and I2 are disjoint open intervals, 
there exist at least one initial condition 7i(a;*) that is in neither sets, and hence the corresponding 
solution 75^ (a; ) satisfies limj-^o 7i (2;) = 0. ■ 

To establish that Ii is non-empty, we show that initial conditions at x = x* sufficiently close to the 
X-axis necessary cross it at some Xj^m < x*. 



Proposition 4.2 For all 70 G [P'(0)x*, 0[, there exists Xmin g]0, x*[ such that the solution o/© only 
exists on X E [xmin, x*] and 7i(xmin) = 0. 

Proof. Pick 72(x) = P'(0)x. Then by (|34l) . we have 

^N(^) = > Vx G [0,x1. 



Applying (pTI) . we get 

7i(x) = P'(0)x+ / P[72](z/)i^[7i,72](z/,x)di/>P'(0)x Vx < x* . (37) 



In particular, 71 cannot cease to exist by diverging to —00 at a finite x < x*. We also see from (|37l) that 
there exists Xq < x* with 71 (xq) > P'(0)xo. Using ^[72] > 0, and 71 (x) > P'(0)x in dlB, we 
find 

li{x) > P'(0)x + (71 (xo) - P'(O)xo) - exp(^ J ji^d, 

> — (p'(0)xo + (7i(xo) - P'(0)xo)e™^"^'W^') . 
Xo V y 

The proof is completed since 7i(xo) — P'(0)xo > 0. ■ 

Before proving that the interval I2 is non-empty, we can establish an additional property of the asymp- 
totically free solution 7*: 

Corollary 4.3 The solution •y^ of Theorem \4. 1 1 satisfies 

7+(x) < 7i(x) < P'(0)x. (38) 

for all X g]0, x*]. 
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Proof. To get the lower bound in (|38l ). we recall that solutions can have at most one local maximum in 
]0, X*] and no local minimum. Since lim7j'(x) = 0, 71 (x) > 7^(x) for all x g]0, x*]. The upper bound 

follows immediately from Proposition I4.2I ■ 

We now show that initial conditions sufficiently close (but above) the nuU-cline 7^(x*) necessarily 
cross it at some x < x*. 



Proposition 4.4 There exist ei ^ 1 sufficiently small such that the solution 0/ (O with 71 (x*) = 
7^(x*)(l — ei) satisfies 71 (x) = •ydx) for some < x < x*. 



Proof. Let < ei < ^ and 71 (x*) = 7^(x*)(l — ei). By continuity of solutions and since 'j^i^*) < 
7i(x*) < 7i(x*)(l — 2ei), there exists < Xi < x* such that 



(39) 



7+(x*) < 7i(x) < 7+(x*)(l — 2ei) Vx g]xi, x*] and 
= iti^*) or 7i(xi) = 'yt{x*){l - 2ei) 
Using these inequalities and hypothesis H2, we get 

d7i . 2ei P(x-) 

-1- > ^ TT^ Vx G Xi,X , 

dx 1 — 2ei7+(x*) 

which, upon integration, gives 

71 (x) < 7+(^^)(l - ^1) + , ^}fl In (-) VxG[xi,xT. (40) 

l-2ei7+(x*) \xj 

Let now X2 be the value at which the r.h.s. of (|40l) attains 7+(x*)(l — 2ei), namely 

/_^+(^^)2(i_2ei)\ ^ 

X2 = X exp — , / , , ] < X . 

2|P(x^)| J 

We now consider the two alternatives X2 < xi and X2 > xi. 

Assume first that X2 < Xi. Then (l40l) shows that 71 (xi) < 7^(x*)(l — 2ei), and thus by definition of 
Xi (see (|39|)). we have 71 (xi) = 7c^(x*) and since by H2, 7^(x) increases as x — 0, there exists an 
X G [xi,x*] with7i(x) = 7^(x). 

Consider now the other possible case, namely X2 > Xi. Since now 7^^ (x*) < 71 (x) < 7^(x*)(l — 2ei) 
for all X G [x2, X*], and we conclude from (|40l) that 

71 (x) - 7e+(x) < itixl - itix) - eati^n + 1 1\ ^^^ff) In ) • (41) 

The proof is completed by noting that for x very close to (but strictly less than) x*, 7^(x*) — 7^(x) is 
negative, while the last two (positive terms) in (|4T]) can be made arbitrarily small by picking ei small 
enough. ■ 
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4.2 Behavior towards x = of non-asymptotically free solutions 



We first consider solutions of ^ corresponding to initial conditions 71 (x*) < 7*(x*). We have the 
following result. 

Proposition 4.5 Let 70 < 7^ (x*). The corresponding solution 71 (x) o/© exists for all x G [0, x*] and 
satisfies 71 (x) = — 1 + C(xln(x)) a* x — 0. 

Proof. Note first that since 70 < 7i(x*), there always exists Xq < x* such that 71 (xq) < 7^(xo). 
Namely, if this does not already hold at x*, the proof of Theorem 14.11 show s that it will eventually hold 
at some smaller value of x. Since solutions can have at most one local maximum and no local minimum 
in [0, X*], we then have 

min(-l,7i(xo)) < ii{x) < 'yti^o) ^ [0>^o] • 
In particular, these solutions exist for all values of x G [0, x*]. We then apply (|301 ) and get 

-1 + Cx+ / Pd2;<7i X < -1 + Cx+ \ ^dz, 

\lt{^^)\jx max(l, |7+(xo)|) A 

where C = i±IiM. By hypotheses HI and H2 (see also dM]), P(x) = P'(0)x + Ofx^) it is then 
straightforward to prove that 71 (x) = — 1 + (9(x(l + ln(x))) as x 0, which completes the proof. ■ 

We now consider solutions of ^ corresponding to initial conditions 71 (x*) in the interval ]7*(x*), 0[. 
We have the following result. 

Proposition 4.6 Let'jo g]7i(x*),0[. The corresponding solution jii^x) of ^ satisfies ■yi^Xuim) = Ofor 
some Xmin £]0, x*[. It can then be continued and enters the first quadrant (becoming double-valued), 
and thus satisfies 71 (xq) > Ofor x > x^am- 

Proof. Let Xmin > be the minimal value such that 71 (x) exists Vx G [xmin,x*]. We claim that 
Xmin > and 7i(xmin) = 0. Namely, since 71 (x*) > 7i (x*), we have by (|3T1) with R[y2] = that 

71 (x) = ^Ux) + iliixn - ll{xn)-^exp( r dz) > 7^(x) (42) 

X \Jx ^7il^j7il^j / 

for all X G [xmin, X*]. Assuming ab absurdum that Xmin = and 71 (x) < for all x G [0, x*] leads to 
a contradiction, for then we would have 71 (x) > 7j'(x) > 7^(x) for all x G [0, x*], and using (|35l) and 
(I42l) . we get 

71 (x) > 7t(x) + (71 (x'^) - il{x'))^^exp(^£ i^^^) 
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which goes to +00 as x — 0. So Xmin > and 7i(xmin) = 0. Akhough ^ is singular at 7i(xmin) = 0, 
these solutions can be continued in the first quadrant by reverting to the so-called 'running coupling' 
formulation of ([5]) (see also (|25l) and [5j). Namely, we introduce a new independent variable t, and 
write X = X{t) and 7i(X(t)) = 71 (t), getting 



dX{t) 



X(t)7i(t) X{to) = X, 



These equations are not singular at 71 = 0, and thus solutions will exist (at least locally around t = to). 
Since P{x) < 0, the solution to these equations will satisfy 71 (t) = 70 > and X{t) = xq > Xmin for 
some finite t > to. ■ 



4.3 Behavior as x ^ 00 

We first show that solutions in the first quadrant are global, and satisfy appropriate estimates as x — >^ 00. 

Proposition 4.7 Let 71 (xq) > 0, and assume P{x) satisfies H3. The corresponding solution 71 (x) 
exists for all x > xq, and satisfies 

X 

< xS'p(xo,x) < 7i(x) < xS'p(xo,x) H 1 

Xo 

for all X > Xq. 



Proof. We first note that 



2 dx V / dx X X ~ 



X 



By integration, we find 



li{x) > xJ + 2 r — = x^p(xo, x) > . (43) 

y ^0 Jxo ^ 

This shows that solutions cannot cease to exist by reaching 71 (x) = at some x > xq. Inserting (l43l) 
into (l30l) gives 

71 (:^) < ""^^ ^ ^'^'"'^^ -1 + x r . dz = x5p(xo, x) + - - 1 , (44) 

Xo Z-^Sp{Xo,Z) Xo 

since 

dS'p(xo,x) P(x) 



dx x^S'p(xo, x) 
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The proof is completed since (l44l) shows that solutions cannot cease to exist by diverging to oo at a 
finite X > Xq either. ■ 

We now turn to the fate of any type of solutions of ^ with 71 (a;*) < as x ^ 00. Our first result is 
that these solutions are global, i.e., they can be extended as a; — > 00. In particular, the asymptotically 
free 7i (x) is global. 



Proposition 4.8 Let 71 (x*) < and assume P(x) satisfies H3. The corresponding solution o/© 
exists for all x > x*, and satisfies 

—xSp{x*,x) < 7i(x) < max [ 7i(x*), sup ) < 

y z&[x*,x] 4 J 

for all X > X*. In particular, ifV{P) < 00, solutions grow at most linearly as x ^ 00. 



Proof. For the lower bound, note first that, as in the proof of Proposition 14. 6 [ we have 

2 dx V / dx X X X ' 

which gives 71 (x) > —xSp{x*, x) upon integration. This shows that solutions cannot diverge to —00 
at a finite x > x*. Now suppose ab absurdum that there exists Xmax < 00 such that 7i(xmax) = 0. By 
hypothesis H1-H3, we have —\ < sup^-g^^,* ^.^^^^^j P(x) < 0, and thus 



7max = max I 7i(x ), sup ^ 

satisfies —\< 7max < 0. Note then that /(71, x) is strictly negative along the 71 = 7max line since 

sup x/(7max, x) < ( 1 + 7max —) < - J ■ 

xe[x*,XmE^^] \ 7max/ 4 

Since 7i(x*) < 7max, this shows that 71 (x) < 7maxforallx G [x^jXmax], contradicting the aZ^afcMrJMm 
assumption. Hence solutions exist globally as x ^ 00. ■ 

Our second result concern the asymptotics of some of these solutions as x 00. Namely, we can 
estimate the growth of solutions that are somewhere less than —1. 

Proposition 4.9 Assume P(x) satisfies H3 and 71 (x*) < 0. If the corresponding solution o/© satis- 
fies 7i(xo) < —Ifor some Xq > x*, then 

X 

—xSp{xo,x) < 7i(x) < — xS'p(xo,x) H 1 < 

Xo 

for all X > Xq. 
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Proof. The lower bound is already contained in Proposition |4]8l The upper bound then follows imme- 
diately from the lower bound and the integral formulation (|30l ). ■ 

The condition 71 (xq) < — 1 is essential in Proposition |4]9] to guarantee that the upper bound is indeed 
negative. We now give possible scenarios that guarantee solutions indeed reach 71 = —1. 

Proposition 4.10 Assume one of the two following statements holds: 

1. — 1 < 7i(x*) < and P{x) satisfies SI or S3, 

2. —1 < 7i(x*) < 7i(x*) and P{x) satisfies S2. 



Then there exists Xq > x* such that the corresponding solution 71 (x) satisfies 71 (xq) = — 1. 



Proof. We consider the alternative 1. first. Note that SI implies S3, so we can use hypothesis S3 only. 
As shown in Proposition 14. 8[ any solution starting at 71 (a;*) < exists for all values of a; > x*. In 
particular, 71 (xc) = 7min < 0. If 7min < —1, the proof is completed. If — 1 < 7inin < 0, we assume 
ab absurdum that 71(0;) > — 1 for all x G [xc, Xd\. Note that there cannot be an x G [xc, x^] such that 
7i(x) = 0, hence — 1 < 7i(x) < for all x G [xc, xj, and we have 

d 2^ o 7i(a:)+7i(3:)' + i 4P(x) + 1 4P(x) + 1 

— (7i(x) ) = 2 ^ > 

dx X 2x 2x 

for all X G [xc, Xd\. Upon integration, we thus find that 



^' 1 + AP{z) 



by hypothesis S3, which is a contradiction. 

Consider then the alternative 2. Under hypothesis S2, we can extend (l35l) to get P(x) < P'(0)x for 
all X G [0, Xc]. Consider now 72(x) = P'(0)x. We have P[72](x) > for all x G [0, xj. Thus, since 
7i(x*) < 7i (x^^) and 7* (x*) < P'(0)x* by Corollary 14. 3 [ we find from (|3TI) that 

7i(x) = P'(0)x- 7i(a;*)-72(x^) i^bi, 72](x^ x) - /p[72](z/) i^[7i, 72](z/, x) dy < P'(0)x 



for all X G [x*,Xc]. The proof is completed since Xc > ~pr^ by hypothesis S2, and thus 71 (xc) < 

P'(0)Xc < -1. ■ 

We conclude this section by showing that 'D{P) < 00 implies that all solutions that are either positive 
or go below 71 = —1 have a finite slope as x ^ 00. 
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Proposition 4.11 Assume V{P) < oo and 71 (xq) < —1 or 71 (xq) > 0, there exists s > such that 

, 7i(x) f -s < if -iiix^) < -1 

hm = < 

X y s>0 if 7i(xo) > 

7/'7i(xo) < — 1, the convergence towards the limit is given by 

-Piz) 





+ s 


poo 


X 




J X 



IfV{P) < 00 and 7i(xo) > 0, then (0) also hold, with —s replaced by s. 



(45) 



Proof. Note that from Proposition 14.71 and 14.91 the hypothesis ^{P) < 00 implies that all solutions 
under consideration here satisfy 



Cl X < |7l(x)| < 02 X 

for some ci, C2 > and all x > xq. The integral formulation (l30l) then gives 

71 (x) l + 7i(a;o) 1 r Piz) 



X 



Xq 



dz 



(46) 



(47) 



from which the proof follows immediately, since the r.h.s. of (|47l) converges by (|46l) and the hypothesis 
V{P) < 00. ■ 



4.4 The confinement solution 

In this section, we consider P{x) satisfying the hypotheses H1,H2 and S4. In particular, recall that we 
assume the existence of x^ > such that P{xr) = —j and P(x) > — ^ for all x > x^ and P(x) tending 
to a finite limit as x ^ 00. Any solution of (|5]) that satisfies 



_ 1 + y/l + 4Poo 

hm 71 (x) = 7oo = ^ , (48) 

x—*oc / 



needs to solve the integral equation obtained by taking the (formal) limit x* — 00 in (|30l) , namely 



Defining 



^^^ ^ k t^{loo + h{xt)f' 700 71 t^{loo + h{xt)r 
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we see that any confinement solution can be written as 7j(x) = ■joo + h{x) where h{x) satisfies h{x) = 
T[h]{x). Consider then B^^ the Banach space obtained by completing the space of C^{[xo, oo), R) 
functions under the norm 



xo = sup +x\f'{x) 

x>xo 



Since lim P{x) — Pqo = lim xP'{x) = 0, \\P — -PoolUo can be made as small as one likes by taking 

X— >oo X— >oo 

Xo > Xr large enough. Standard arguments then show that T is a contraction in a ball of positive 
radius p < 1 centered at in B, which shows there exists a unique h E B solving h = T[h]. Since 
h{x) is regular, 7^(2;) = + h{x) solves ([5]) for all x > xq. We now remark that 'Ji{x) will satisfy 
(|49] ) as long as it exists when x decreases below xq. However, it can only cease to exist if it satisfies 
7i(a;min) = for some Xmin > 0, for the r.h.s. of © is negative for large negative 71. Assuming it can 
be continued up to x = 0, we can conclude from Theorem (I2.2|) that either 7i(0) = or 7^(0) = — 1. 
Finally, if P'(x) > for all x > Xmax, then the lower nuUcline 



_ ^l + AP{x) + l 
7i [x) = ^ 

is decreasing towards 700- If there was an xq > Xma.x such that ^i{xo) = 7f (xq), then 7i(a;) would 
enter a region of strictly positive derivatives w.r.t. x, and hence we would get 7^(0;) > 7f (xq) > 'joo, 
contradicting (|48l) . Similarly, if there was an xq > Xmax such that 7i (xq) = 7ooj then 7^ (x) would enter 
a region of strictly negative derivative w.r.t. x, and hence would satisfy 7j(x) < 700 for all x > xq, 
contradicting (|48|) again. We thus find that under the monotonicity assumption on P(x), we have 



1 + VI + 4P00 , , 1 + v/l + 4P(x) 
7oo = ^ < < 



2 ' ' ' 2 

for all (finite) x > Xmax- This concludes the proof of Theorem 12. 51 



5 Post Scriptum: QED[5] revisited 



In a previous publication (see [|5]|), we have considered the initial value problem 

d7i(a;) .^ . _ 7i(a;) +7i(a;)^ - P(a;) , . 

= fs{li{x),x) = , 7i(xo) = 7o > . (50) 

dx sx7i(x) 

with < Xo < 1 fixed and P a function on [0, 00), positive on (0, 00), with P(0) = and P'(0) 7^ 0. 
Defining the following (possibly infinite) quantity 

we showed that 'D(P) was intimately linked to the behavior/existence as x ^ 00 of solutions starting 
with 7i(xo) > 0. Namely, if ^^{P) < 00, we showed that there was a smallest (non-zero) value 7i (xq) 
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separating global solutions (with 71 (xq) > Ji{xo)) from solutions that cannot be continued to x = 00. 
We also showed that despite the singular nature of (|50l ) as x — > and/or 71 — > 0, all solutions of (|50l ) 
could be continued for all x E [0, Xq], and would approach (0, 0) while satisfying for all x G [0, Xq] the 
bound 



Cb X if s < 1 

7i(a;) < { Cb X \ ln(x)| if s = 1 , 
Cb x^l' if s > 1 



(51) 



for some constant C^ = Cb{s, 70). In Lemma [5^ below. we will improve the above to get linear bounds 
in all cases (i.e. for s > 1 as well). In the mean time, we define, for any (finite) integer p > 2, the 
truncation of the (divergent) series solution: 



72 



,p(x) = P'(0)x + J]a„ x". 



n=2 



The coefficients {an}n=2...p can be found recursively by imposing |i?[72,p](x)| < Cx^ ^ as x 
where the remainder map R is defined by 



0, 



dx 



SX72(x) 



for 72(x) any function on [0,xo]. Note that by choosing xq sufficiently small, we can ensure that 
72,p(x) > for all x G [0,xo]. Also, since -P'(O) 7^ and 72,^ is continuous, there exists a constant 
C > such that 72 (x) < Cx for all x G [0, xq]. 

Finally, we also note that if 71 solves (|50l) and 72 (x) is any positive function on [0,Xo], then for all 
X G [0, Xq], we have 

7i (a;) - 72(3;) =(71(3^0) - 72(3:0)) ^[71, 72](a;o, a;) + / R[y2]{y) K[yi,-f2]{y,x) dy , (52) 

where 

^[7i,72](a;o,a;) = exp( 



1 



+ 



P(z) 



xo SZ SZ-ii{z)-f2{.z) 



dz 



exp 



xo 



P{z) 



SZ^l{z)-l2{.Z 



-dz 



We are now in position to show that all solutions to (1501) agree to all orders in perturbation theory with 
the series solution. 



Theorem 5.1 Let 71 and 72 be two solutions o/fl50l). Fix p > 2 and let 72 „ as above. Then there exist 



constants C and Cp such that 



li{x) - 72,p(x) I < CpX^ , 
liix) - 72 (x) < 71 (xo) - 72(a;o) 



X 

— J exp 

Xo 



Vxo X 



(53) 
(54) 



where (122) holds for all x G [0, xo/2] and ({54h for all x G [0, xq]. 
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In other words, for any solution of (|50l) . any derivative (of finite order) converges as x — > to the 
corresponding derivative of 72,^. Hence a (truncated) power series expansion at x = of any solution 
agrees to any (finite) order with the truncated divergent series of the same order. Also, the difference 
between any two solutions of (l50l) decays faster than e"*^/^' as x ^ 0. 

Proof. In Lemma [5^ below, we prove that there exist constants ci and C2 such that 71 (x) < cix and 
72(x) < C2X for all x G [0, Xq] if 71 and 72 are solutions of (|50l ). In particular, there exist a constant 
C > such that 

K\r.M(y,.) < (^) ' exp(- £|d.) = (£) '• exp(c(i - i)) (55) 
for all < X < ?/ < Xq. Since -R[72] (x) = if 72 is a solution of (l50l) . we have 

71 (a;) -72(3;) =(71(2;) -72(a;))/'sr[7i,72](xo,x) , 
from which we get (l54l) immediately. 

On the other hand, 72,p also satisfies 72,^(2;) < Cx for all x G [0, xq], and so K[yi, 72,p] also satisfies 

i^[7i,72.](y,a:)<(^)%xp(c7(i-i)) 

for all < X < ?/ < Xq. In particular, for fixed xq, K[yi, 72,p](xo, x) decays faster than any power law 
as X — s> 0. Assume now that x < Xo/2, and note that 

I -72,p(a;) I < I 71(3^0) -72,p(a;o) | ^^[71, 72,p](a;o, a;) + C y^'^ K['ji,'j2,p]{y, x) dy . 



Splitting the integral over [x, xq] into [x, 2x] and [2x, xq], and using that K[yi, 72,p] {y, x) is a decreasing 
function of ?/, we find 



2x r2x 



y^-'K[^i,l2,p]{y,x)dy< / yP-iir[7i,72,J(x,x)di/ = x^(^) 

J X 

yP-^ if[7i,72,p](i/,x) dy < ^[71, 72,,](2x, x) / y^-^ dy = K[^,,^,,^]{2x,x){^-^) . 



2x J2x 



The proof of (l53l) is completed by noting that by (l55l) . K[yi, 72,p](2x, x) and K[yi, 72,p](xo, x) decay 
faster than any power law as x ^ 0. ■ 



Lemma 5.2 For any 70 > 0, the solution o/(|50l) exists for all x G [0, xq] and 

7i(x) <Cx Vx G [0,xo] 

for some C = C(xo,7o) > 0. 
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Proof. We first note for future reference that is continuous for all x G [0, Xq], and that there exists 
constants C± > such that C_x < P{x) < C+x for all x E [0, Xq]. By (|5T| ). we only have to consider 
s > 1. We first choose a > such that 



a > 



max(2C+, Cb| In(xo)l) if s = l 

max(2C+,^^;^,Q4-i(s-iri) if s > 1 



We then note that 



fs{a x,x)-a = — - ^ + a(i _ 1) > ^ _ + a(l _ 1) > ^ + a(i _ 1) , 

sx sx^OL sx sxa 2sx 

from which it follows that 

( xo if s = 1 

fs{ax,x)>a y X E [0,X{s,a)] where X(s, = ^ if s > 1 ' ^^^^ 



4q:(s-1) 



We then note that 



CbX{l,a)\\n{X{l,a))\ = C;, Xo| ln(a;o)| < axo = 

aXis,ay/^ = < 4(i^ = if s>l. ^^'^ 

In other words, (|5T| ) and (|57] ) imply that the solution of (|50l ) satisfies 7i(X(s, a)) < aX(s, a). Since 
X(s, a) < Xq, this means that as x decreases, solutions enter the triangle 

'^s,a = { (x, 7i) I X G [0, X{s, a)] and < 71 < a x } 

through its right boundary (i.e. at x = X(s, a)). 

Solutions that enter „ at x = X{s,a) cannot satisfy 71 (x*) = ax* at some x* < X(s, a), for 
by (l56l) . we would have 71 (x) > ax for all x G (x* , X{s, a)], a contradiction with 7i(X(s, a)) < 
aX(s, a). Hence 71 (x) < ax for all x G [0, X(s, a)], and the proof is completed by noting that the 
bound (|5T|) is stronger than 71 (x) < ax for x G (X(s, a), Xq]. ■ 
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